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Most helix-coil transition theories can be characterized by three parameters: energetic, describing 
the (free) energy cost of forming a helical state in one repeating unit; entropic, accounting for the 
decrease of entropy due to the helical state formation; and geometric, indicating how many repeating 
units are affected by the formation of one helical state. Depending on their effect on the helix-coil 
transition, solvents or co-solutes can be classified with respect to their action on these parameters. 
Solvent interactions that alter the entropic cost of helix formation by their osmotic action can 
affect both the stability (transition temperature) and the cooperativity (transition interval) of the 
helix-coil transition. A consistent inclusion of osmotic pressure effects in a description of helix-coil 
transition, for poly(L-glutamic acid) in solution with polyethylene glycol, can offer an explanation 
of the experimentally observed linear dependence of transition temperature on osmotic pressure as 
well as the concurrent changes in the cooperativity of the transition. 
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The helix-coil transition is central to many processes 
in living matter [HQ- To mimic the interactions and 
structures found in nature, in vitro experiments related 
to biopolymers are usually performed in solutions of dif- 
ferent composition. Naturally, a theoretical description 
of solvent composition effects is a necessary component 
of any helix-coil transition theory. There are several ap- 
proaches that offer such a description within the standard 
Zimm-Bragg model @, i|. Farago and rincus [oj| pro- 
posed a classification of solvents based on their action on 
the parameters within this model. They have shown that 
solvents that promote the relative stability of the helical 
state compared to the coil rescale the parameter s, alter- 
ing the melting temperature, while solvents that affect 
the helix-coil intcrfacial free energy modify a, changing 
the melting interval. Recent circular dichroism measure- 
ments of the helix-coil transition of poly(L-glutamic acid) 
(PLGA) in an aqueous solution of co-solvent polyethylene 



glycol (PEG) |6| show that besides stabilizing helices, 
which agrees with stability-altering solvent description 
of Farago and Pincus, the increase of PEG concentration 
also couples with the decrease of cooperativity, which is 
harder to explain within the Zimm-Bragg approach or its 
modifications Q . The main action of PEG is to modify 
the solvent-accessible volume within the polypeptide (sec 
Fig- [J), an d can thus fundamentally change the polypep- 
tide conformational space and the associated conforma- 
tional entropy. It is this coupling between the osmotic ac- 
tion of the PEG and the modifications on the polypeptide 
conformational space that will be elaborated in what fol- 
lows. Here we propose a consistent theoretical framework 
in which a single osmotic action of PEG solution, encoded 
in the osmotic pressure variation of the free energy dif- 
ference between the helical and coil repeating units, en- 
genders changes both in the melting temperature and in 
the melting interval, thus correctly describing the exper- 
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FIG. 1. Two fragments of 15 repeat units (1-15 of 2JU4 and 
16-30 of 1BBA) from Protein Data Base structures are shown 
to illustrate the change in solvent accesible volume AV be- 
tween disordered and a-helical conformations. 



imcntal data. 

A minimal formulation of the Zimm-Bragg model im- 
plies a characteristic equation of the form (1 — A)(s — 
A) — sa = 0, where the two phenomenological param- 
eters are defined as s = exp( AH ^F ), and a is the 
temperature-independent entropic cost of creating a he- 
lix domain within coil regions |J|. The difference of en- 
thalpies AH = Hheiix — H CO ii accounts for the relative 
energetic gain of helix formation, and AS = Sheiix — S co u 
describes the entropic cost of helix formation, with T 
the absolute temperature, measured in units of ks- The 
helix-coil transition occurs at the point where the en- 
ergetic gain is compensated by the entropic loss, i.e., 
s = cxp(AiJ/T)/ cxp(A5) ~ 1, resulting in the melting 
temperature T m ~ AH/ AS. The interval of transition 
can be estimated as AT ~ s/a. Within the Zimm-Bragg 
formulation, the solvent effects that alter the parameter s 
would lead to changes in the stability of the helix, mod- 
ifying its transition temperature T m , while the solvent 
affecting a would alter only the cooperativity of tran- 
sition AT @. In this context, s and a are introduced 
as completely independent parameters corresponding to 
two different mechanisms of solvent action. Any solvent 
effects that simultaneously affect both transition stabil- 
ity and cooperativity would then require more than a 
single mechanism of action, which is hard to imagine for 
such simple solute as PEG. Even within this Zimm-Bragg 
parametrization, one might suspect that the relative en- 
tropic cost of helix formation AS and the entropic cost of 
domain formation a could be in principle related, since 
they both include the entropy difference between the re- 
peating units in the helical and coil conformations. 



The helical structure of biopolymers is mainly stabi- 
lized by intermolccular hydrogen bonding between re- 
peating units, the presence of hydrogen bonds being a 
prerequisite for helix formation. A statistical descrip- 
tion of the solvent effects on the helix-coil transition re- 
quires at least three parameters: an energetic parameter, 
W = V + 1 = cxp(C/ /T), where U is the energy of a hy- 
drogen bond; an entropic parameter, Q, that stands for 
the ratio between the number of all accessible states ver- 
sus the number of states available for the repeating unit 
in a helical conformation; and a geometric parameter, A, 
that describes the sequential geometry of hydrogen bond 
formation. We first describe the helix-coil transition in 



35 





W 

Q 


a) AE > V_ 
Q<Q ^ 




1.1 1.2 


1.3 


b) AE = \ 
Q = Q \ 


W 

Q " 




0.9 1 


l.i 


c) AE < A 
Q>Q 


w 

Q 



0.8 0.9 1 

t=T/T 

FIG. 2. Temperature dependence of two largest eigenvalues at 
different AE = E h - E c ; Q = 60, A = 3, q = 3. T = U/ In Q 
is the melting temperature of solvent-free model. 

vacuo, i.e., without any solvent. The Hamiltonian of a 
solvent-free model [1, [9[ can be written in the form 

N 

-0H O (H}) = JE^ A) 

as considered in @. Here /3 = T -1 , N is the number of 
repeating units, and J = U /T is the temperature-reduced 

energy of hydrogen bonding, <Sj A) = 11^=0 s ilj+k, 1), 
where 5(x, 1) stands for the Kroncckcr symbol and 7; = 
1,...,Q. The spin variable 7 describes the state of each 
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repeating unit by assigning to it one of the Q possible 
conformations: conformation 1 corresponds to the heli- 
cal, and the remaining Q — 1 conformations to the coil 
state. The partition function 



Q N 

z (v,q)= ^nl 1 



VS. 



(A) 



{7*=l}i=l 



(2) 



can be evaluated by applying the transfer matrix method, 
corresponding to Eq. Q] (see Ref. resulting in the 
characteristic equation 



A A ~ X (A - W)(X -Q) = (W- 1)(Q - 1). 



(3) 



In the thermodynamic limit, the problem simplifies con- 
siderably. It is enough to study the two largest eigenval- 
ues obtained from Eq. [31 Numerically they are closest at 
the point where asymptotes W(T) and Q cross (Fig. 0b), 
a point defining the transition temperature To. This is 
in accord with general physical considerations: transi- 
tion occurs when entropy and energy compensate each 
other, at To = U/hiQ. The minimal distance between 
the two eigenvalues can be estimated as Q 1_A (see Q) 
and is related to the helix-coil transition interval, a co- 
operativity measure. To quantify the transition interval, 
it is appropriate and informative to introduce the spatial 
correlation length £ = In -1 (xjj! where Ai and A2 are 
the first and second leading eigenvalues of the character- 
istic equation. The temperature dependent £(T) has a 
maximum at the transition point. Its value is related to 
the transition interval as AT ~ £max ~ Q^~^ 0- 

Up until this point we have not taken into consider- 
ation any specific co-solute or solvent effects. Solvents 
able to promote solvent-polymer hydrogen bonds have 
been analyzed in Refs. @ and 11, Here, however, we con- 
sider another vast group of co-solutes which do not affect 
hydrogen bonding directly, but do modify the polypep- 
tide conformations by changing the chemical potential 
of the solvent or the osmotic pressure of the solution. 
A classical example of such a co-solute is PEG, which 
can act as an osmoticant and as a depletion agent fl3| . 
Because of their size, PEG molecules are depleted near 
the polypeptide chain, exerting an osmotic pressure that 
changes the energetic cost of certain conformations vs. 
others. These effects of PEG are well documented and 
have been explored extensively [l6|. The osmotic pres- 
sure of the solution depends only on the concentration of 
PEG, provided that all the other components can equally 
access the helix and the coil state of the polypeptide. It 
is a known function of its concentration that has been 
investigated in detail [l7j |. 

Our model of solvent is based on the following as- 
sumptions: i) solvent molecules can interact with (affect) 
both helical and coil repeating units of the polypeptide; 
ii) binding of solvent molecules changes the free energy 
of a repeating unit, depending on the conformation of 



the repeating unit (Eh for helix, and E c for coil); iii) 
polypcptide-solvent interaction depends on the orienta- 
tion of solvent molecules around the repeating unit and 
the number of solvent orientations q > 2 account for the 
solvent entropy; iiii) a spin variable fii € [l,g] describes 
the state (orientation) of a solvent molecule and its value 
1 corresponds to solvent binding. 

While the various contributions to Eh, E c are difficult 
to disentangle, the overall difference AE = Eh — E c can 
be analyzed explicitly. The (free) energy difference Ai? 
describes the effect of the co-solute, and the larger this 
difference, the stronger the stabilization of the helical 
state vs. the coil state. Since the main action of PEG is 
depletion-induced osmotic pressure, this free energy dif- 
ference or cquivalcntly the corresponding osmotic work 
required to drive the chain through the helix-coil transi- 
tion can be written as 14 



\AE\ 



n n 



AV. 



(4) 



Here n osm its the osmotic pressure of PEG solution, 
and Ay is the volume of water that must be exchanged 
with the bulk when the polypeptide chain goes through 
the transition [HI (see Fig. Q]). This means that the 
helix-coil transition responds in a way analogous to a 
semi-permeable membrane, excluding solvent from cer- 
tain portions of the polypeptide chain. The sign of AE 
depends specifically on the details of the osmotic action 
of PEG, stabilizing or destabilizing the helix vs. the coil 
state. 

We now add these solvent-mediated changes to the free 
energy of a repeating unit in the original in vacuo Hamil- 
tonian Hq of the helix-coil transition [1 - 11 1 . This results 



-0H ({j h fit}) = -pH ({ 7i }) - (3H solv ({ 7i }, {m}) 



N 

= E H A) + J c (1 - s (m, 1) + hs^s fa, 1)) , 

i=l 

(5) 

where I Ct h = E Ct h/T. The second and third terms on the 
rhs of Eq. [5] describe solvent interactions with repeating 
units in the coil and helical conformations, respectively. 
The partition function is then 

Q Q 

Z(V,I h ,I c ,Q)= E expi-pHd^m})). (6) 

{ T i=l}{Mi=l} 

Each of the solvent degrees of freedom, is coupled 
with a single 7, and the lattice expansion for Eq.|6]is of a 
" decorated lattice" type 12( . This allows one to sum out 
the solvent degrees of freedom. After a re-arrangement 
that makes use of the Kronecker delta properties (for 
details see [13), Eq. [S] becomes 



Z(V,I h ,I c ,Q) = (q + e Ih ) N Z (V,Q), 



(7) 
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where 

When Eq. [7] is compared with Eq. [21 it becomes obvious 
that the two expressions are identical up to an irrelevant 
multiplicative factor. Thus we can conclude that the only 
effect of the osmotic action of PEG is the transformation 
of the entropic parameter Q — > Q. Therefore, the same 
characteristic equation, Eq. [3l but with a redefined Q re- 
mains valid: X A - 1 (X-W){X-Q) = (W-1)(Q-1). Since 
Q is now temperature dependent, the W ~ Q condition 
occurs at temperatures different from the melting tem- 
perature of the original model To- In fact, two regimes 
of temperature dependence of Q are possible, as shown 
in Figs. [21 a) and c), depending on the sign of AE. As- 
suming first that the PEG osmotic pressure stabilizes the 
helix, AE < 0, the helix-coil transition then takes place 
at higher temperatures (Fig. [21 a)). However, this is not 
the only effect of the PEG osmotic action: the distance 
between eigenvalues also becomes larger and the transi- 
tion interval thus decreases. For AE > 0, the situation is 
reversed: asymptotes cross at lower temperatures (desta- 
bilisation) and the transition interval increases (Fig. [21 
c)). What is interesting is that even after the rescaling 
of parameter Q — > Q, Eq. the point of closest ap- 
proach between the two largest eigenvalues can still be 
estimated from the crossing between W and Q curves. 

To study the stability and cooperativity of this gen- 
eralized model, we analyze the temperature dependence 
of the correlation length from the computed eigenvalues. 
This allows us to calculate the transition point (T m ) and 
the cooperativity (£ max ) as a function of AE (Fig. [3}. As 
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FIG. 3. The cooperatively of the transition quantified by £ m ax 
as a function of AE ~ n osm . Inset shows the concurrent 
dependence of the melting temperature t m = T m /To. 

shown in Fig. [3] (see inset), the helix-coil melting temper- 
ature grows linearly with increased AE, indicating the 
increased stability of the system. This is consistent with 
the experimentally observed linear increase of the tran- 
sition temperature with osmotic pressure as reported by 



Stanley and Strey @, @] ■ On the other hand, the coopera- 
tivity measure £, ma ,x decreases at the same time (Fig. [3]), 
as has been observed by Koutsioubas et al. [6(. This 
two-fold action of the PEG cannot be captured by the 
Zimm-Bragg model @, @] , in which the stability and the 
cooperatively of the transition are described by two in- 
dependent parameters. 

Why can't the Zimm-Bragg based approach describe a 
decrease of cooperativity concurrently with an increase in 
stability? As shown in Q, the Zimm-Bragg parameters 
can be recast in terms of our parameters as s = W/Q 
and a = Q 1_A ~ £~ 2 . The entropic parameter Q is 
present in both s and a. Obviously its changes will alter 
both the stability and the cooperativity of the system. 
The assumption of independence of s and cr, an inher- 
ent property of Zimm-Bragg and related models, leads 
to a description of the effects of PEG [H that cannot eas- 
ily be reconciled with experiments Ref. @. We showed 
above that the resolution of this discrepancy is impossible 
within the Zimm-Bragg model and that an alternative, 
more detailed microscopic description of the PEG action, 
is necessary, which in its turn leads to a picture consistent 
with experimental results. 
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